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Measurement-induced nonlocality is a measure of nonlocalty introduced by Luo and Fu [Phys. 
Rev. bett 106, 120401 (2011)]. In this paper, we study the problem of evaluation of Measurement- 
induced nonlocality (MIN) for an arbitrary mxn dimensional bipartite density matrix p for the case 
where one of its reduced density matrix, p a , is degenerate (the nondegenerate case was explained in 
the preceding reference). Suppose that, in general, p a has d degenerate subspaces with dimension 
mi (mi < m, i = 1, 2, d). We show that according to the degeneracy of p a , if we expand p in a 
suitable basis, the evaluation of MIN for an m x n dimensional state p, is degraded to finding the 
MIN in the mi x n dimensional subspaces of state p. This method can reduce the calculations in 
the evaluation of MIN. Moreover, for an arbitrary mxn state p for which mi < 2, our method leads 
to the exact value of the MIN. Also, we obtain an upper bound for MIN which can improve the 
ones introduced in the above mentioned reference. In the final, we explain the evaluation of MIN 
for 3 x n dimensional states in details. 

PACS numbers: 03.65.Ud, 03.67.Mn 



INTRODUCTION 

Quantum mechanics is a nonlocal theory The princi- 
ple of locality states that the properties of one particle 
can not be affected by another particle that is sufficiently 
far away. Nonlocality in quantum mechanics, at least, 
has two different aspects On of these aspects arises 
in the Aharonov-Bohm effect. The Aharonov-Bohm ef- 
fect is nonlocal in the sense that the electromagnetic field 
influences an electron in a region where the field is zero. 
The other aspect of nonlocality in quantum mechanics 
appears when one performs local measurements on spa- 
tially separated systems. These nonlocal effects result 
from the fact that the local measurements can alter the 
overall state of a multipartite quantum system. Quantum 
nonlocality usually refers to this aspect of nonlocality and 
is studied often in the context of Bell inequalities. Bell 
inequalities are mathematical relations setting conditions 
on the results of measurements made on separated sys- 
tems. These inequalities are satisfied by any local hidden 
variable theory, but they may be violated by quantum 
mechanics. This is the very feature of quantum mechan- 
ics that is usually mentioned as quantum nonlocality. 

Recently, S. Luo and S. Fu in Refs. 0, Q studied the 
latter aspect of nonlocality based on an approach differ- 
ent from the violation of Bell inequalities. They have 
used this idea that, in general, measurement in quan- 
tum mechanics causes disturbance. If one performs lo- 
cal measurements that do not disturb the states of the 
subsystems, then any disturbance in the system's overall 
state, can be attributed to the genuine nonlocal features 
of the system. The disturbance caused by local measure- 
ments leaving the states of the subsystems invariant can 
be quantified by distance between the overall pre and 
post-measurement states of the system. Based on this 



idea, they defined the measurement-induced nonlocality 
(MIN) as the maximum distance between the bipartite 
state p and its post-measurement state, where the max- 
imum is taken over all the von Neumann local measure- 
ments which do not disturb the local state p a = tr^p 0] . 

The MIN is a manifestation of quantum nonlocality be- 
sides the violation of Bell inequalities. It is different from 
entanglement and discord, although for any pure state it 
coincide with the geometric measure of quantum discord 
and the square of concurrence @, 0] ■ MIN can provide a 
novel classification scheme for bipartite states, and it can 
also be a quantum resource in quantum information pro- 
cessing, although there is no operational interpretation 
for it. 

The analytical formulas of MIN for any state with non- 
degenerate local state p a — tr\,p, arbitrary dimensional 
bipartite pure states and 2 x n dimensional mixed state 
was obtained in Refs. 0, Q. In addition, in Ref. 0], 
a tight upper bound on the MIN of an arbitrary mxn 
dimensional state is derived. The necessary and sufficient 
conditions for a state to have nullity of MIN were studied 
in Ref. In this paper, we study the evaluation of MIN 
for an arbitrary m x n dimensional bipartite state. We 
introduce a method that can reduce the calculations in 
the evaluation of MIN. Using this method, we obtain 
an upper bound for MIN which can improve the ones 
introduced in Ref. 

The paper is organized as follows. In the next section 
(Sec. II), we review the results of Ref. Our main 
results are given in Sec. III. In Sec. IV, 3 x n dimensional 
states are studied. Finally, we give some conclusions in 
Sec. V. 
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MEASUREMENT-INDUCED NONLOCALITY 

Consider a bipartite state p defined in Hilbert space 
H a <g> H h where H a (H b ) is the Hilbert space of part a 
(b) with dimension to (n). Suppose this state is shared 
between two distant parties A and B. Assume that 
party A performs a non-selective local von Neumann 
measurement on his (her) part, then the state of system, 
p, changes to n Q (p) = Efc( n £ ® lfc )p( n £ ® ! fc ) whc rc 
n a = {n^} is a set of orthogonal, one dimensional pro- 
jectors summing to the identity. Then, the measurement- 
induced nonlocality is defined as 

N(p)=xa^\\p-U a (p)f. (1) 

where the maximum is taken over all the von Neumann 
measurements II a = {H k } which do not disturb p a lo- 
cally, that is, Efc n fcP an fe = P a and INI 2 is thc Hilbert 
Schmidt distance defined as ||X|| 2 = trX^X. The maxi- 
mum in definition of MIN is used to capture all the nonlo- 
cal effects that can be induced (indicated) by local mea- 
surements. 

When p a is non-degenerate with spectral decomposi- 
tion p a = ^\ Afe|fc)(fc|, then the only von Neumann mea- 
surement that does not disturb p a is {IIJJ = so in 
this case, the maximum in Eq. (TTJ) is not necessary. In 
Refs. 0, 3, some other basic properties of the MIN, are 
listed. 

Any to x n bipartite state p can be expanded in 
terms of {Xo = I/y/m,Xi : i = l,...,m 2 — 1} and 
{Y = I/ y /fi,Y J :j = l,...,n 2 -l}a8 




where {Xi : i = 1, to 2 — 1} and {Yj : j = 1, n 2 — 1} 
arc traceless Hcrmitian operators satisfying the condi- 
tions trXiXii = Sa' and trYjYy = Sjj> respectively. It 
was shown that 0] MIN can be expressed as: 

N(p) = trTT 1 - nun trAT^A 1 . (3) 

where T = (tij) is an (m 2 — 1) x (n 2 — 1) dimensional 
matrix and A = (aki) is an m x (to 2 — 1) dimensional 
matrix with aki = t r (J^%Xi). Here {n^ = is any 

von Neumann measurement leaving p a invariant. In Eq. 
([3]) "t" denotes transpose of matrices. 

Denning a kQ = tr(\k)(k\X ) = l/Vm, then {a kl : i = 
0, 1, . . . , to 2 — 1} are the coefficients for expansion of the 
operator |fc)(fc| in terms of {Xi : i = 0, ...,m 2 — 1} and 
thus 

m 2 -l 

£ ak l a k , l =tr\k){k\k'){k'\=5 kk ,. (4) 

i=0 



where k,k = 1,2, ...,m. The author in Ref. consid- 
ering only this latter constraint, obtained the following 
upper bound for MIN: 

m 2 — m 

N(p)< J2 A - (5) 

i=l 

where the {A, : i = 1, ...,to 2 — 1} are the eigenvalues of 
the matrix TT , arranged in decreasing order. In addi- 
tion, in Ref. Q , the analytical formulas of MIN for any 
dimensional pure state and 2 (gin dimensional mixed state 
was obtained. 

In the next section, we study the problem of evalua- 
tion of MIN for an arbitrary to x n dimensional bipartite 
density matrix p where its reduced density matrix, p a , 
has d degenerate subspaces with dimension rrii. 

MAIN RESULTS 

In order to evaluate the MIN, at first, we focus on 
the constraints in the optimization problem in Eq. ((T|). 
In this equation, the maximum is taken over von Neu- 
mann measurements leaving p a invariant. A von Neu- 
mann measurement is defined via the set n° = {n^} 
where the set H a = {H%} is a set of orthogonal, one- 
dimensional projection operators summing to identity i. 
e. H£H« = 5 W U%, trill = 1 and £™ x U% = I. These 
relations set some conditions on {aki}- If we expand 
n£ in terms of {Xi : i = 0,1,..., m 2 - 1} i. e. n£ = 

m 2 — 1 

2i=o a kiXi, then from the relation n^nj!, = Skk'^-t 
we deduce: 

m 2 -l m 2 -l 

/~] akiak'i'XiXi' = Skk' 2^ aki -^ i ■ (^) 

i,i'=0 i=0 

In addition, from Y^^—i ^fc = I we have ^2 k a^Xi = 
I = ^/mX®. So from the orthonormality of {Xi} we 
obtain 

m 

^a fc! ;=0, i = l,2,...,m 2 -l. (7) 

fc=i 

Beside the previous constraints, the invariance of p a un- 
der von Neumann measurements dictates the relation 
P a = Efe n fe/° a IIfc- This relation is fulfilled if and only 
if n a = {n^} is an spectral projections of p a , that is 
n a = {H c k = \k)(k\} where \k) are eigenvectors of p a . 
If p a is non-degenerate, then the {\k}} is unique and so 
the maximum in Eq. ([I]) is not necessary. On the other 
hand, when p a is degenerate, any linear combination of 
eigenvectors corresponding to the same eigenvalue, is also 
an eigenvector of p a . But still, for every non-degenerate 
eigenvalues of p a there is a unique eigenvector. So, for an 
arbitrary p, taking maximum over non-degenerate sub- 
spaces of p a in Eq. ([T]) is not necessary and only degen- 
erate ones contribute to the maximum in this equation. 
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If the constraints in Eqs. ([5]), ([7]) can be divided into con- 
straints for each degenerate subspacc of p a , independent 
of the other subspaces of p a , then the maximum in Eq. 
((T|) is degraded to the sum of the maximums for every 
degenerate subspaces of p a separately. This can be done 
by a suitable construction of {Xi} from the eigenvectors 
of p a as we do now. 

We expand an arbitrary m x n dimensional state p in 
the form 

m 2 n 2 -l 
i=l j=Q 

where the set {Yj} are the same as the one used in Eq. @ 
and the set {X^} is constructed from the eigenvectors of 
p a as follows. Let us first introduce the notation we use in 
the next relations. We denote the spectral decomposition 

Of p a by p a = Yf r=1 e r J2T=1 \krs}(krs\+T,t=l e 'r\ k r)(k r \ 

where {e r : r = l,...,d} ({e' r : r = l,...,d'}) are the 
degenerate (non-degenerate) eigenvalues of p a . In this 
relation we assume that p a has d + d!(l < d + d 1 < m) 
different eigenvalues and each eigenvalue e r is repeated 
ra r times (2 < m r < m,d' + Er=i m r = m )- Also, the 
degenerate eigenvectors of p a are denoted by \k rs ) where 
the first index specifies the corresponding eigenvalue and 
the second one is related to the degeneracy of the eigen- 
value specified by the index "r". The non-degenerate 
eigenvectors of p a are represented as \k r ). 

At first, we recall that the set {|fc rs )} is not unique 
and any linear combination of \k rs ) with the same eigen- 
value is also an eigenvector of p a . For the construction 
of {X-}, we choose one of these sets which we denote as 
{\k' rs }}. We now construct from {\k' rs )(k' rs ,\ : s,s' = 
l,...,m r } the operators {X- : i = B r -i + 2,...,£> r }, 
where B r = m 2 + ... + m 2 and Bo = 0, in such a 
way that they form a set of traceless Hcrmitian op- 
erators satisfying the condition trX^Xj = 6{j and we 
set X'l+Br-! = ^T,T=i\Ks)(Ks\- Moreover, we put 
{X' Bd+r = \k r )(k r \} and the set {X[ : i = D,...,m 2 } 
(Bd + d' + 1 = D) arc constructed from the mixture of 
\k' rs )(k r > | such that the relations X[ = X i , trX[X'^ = 5ij, 
H\K s )(Ks\ X i) = and tr(\k r ){k r \Xr) = are satisfied. 

Now, we evaluate MIN using the expansion of p in the 
new basis. As stated before, the IP = {II£} used in Eq. 
(p} must be an spectral projections of p a i. e. in our 
notation n a = {Ilj^ = |fcp»>(* r .|} + {ng p = \k r )(k r \} 0. 
After performing von Neumann measurement 
IF leaving p a invariant, the state p changes to 

n» = EtiEr^nL^Mn^iHEtiOi^ ® 

l)p(H% <8> 1). Using the Eq. © we can write 

n» r = Eti E^Ti E£ Eji 1 c ij u ts x i u k r s ® 

Yj + Eli E™i E"=o 1 dPtXPt ® Y r From thc 
construction of {X[} we have {IK, X'^% = : i ^ 

F r , B r } and {Eg X<IL% = : * j?B A + 1, B d + d'} 



where F r = B r _\ + 1. Thus, we conclude that 

n» = EtiErJiEf^ES'^n^^'n^ 

Yj + EKi E^o 1 ® Therefore p - IP(p) = 

EtiE^ESo 1 ^-^ E^n^n^) ® 

+ E"=i5 Ej=o Cij^i' ® ^j- On thc other hand, 
from Eq. ([8]) and the invariance of p a under the von 

2 

Neumann measurement we get p a = Ei=i c i§X[sJn 

and P a = EtiEr=i n L/>° n L + X)£=i ng^ng, 

2 

respectively. Consequently, we have Ei=i c io^' = 

Ed Tr^m r Tr^B r no Y'Tta , ir^B d +d' v , 

r=ll^ S =ll^i=F r C i0^k r3 A i ll k r3 + l^i=B d +l C iO A i 

and thus EtiEf= d F r Ciopq - ££x n^ngj + 

2 

E™ _d c i0-^i = 0. Therefore, we conclude that in the 
latter relation for p — IP (p), the term j — vanishes. 
Now, we define K rs i := <rll^ X^, then 

n 2 -l m 2 

||p-n Q ( P )|| 2 = E(E4 

j=l i=D 
d B r B r m r 

+EE (4 -EE ^-fifrri' en *,j )) (9) 

r— 1 i=F r i' = F r s— 1 

The Eq. ^ and Eq. ([7]) are changed to 

B T B r 

^ ' K rs iK rs ni XiXi' = <5 SS ' S ' K rs iXi, (10) 

i,i'=B r _i+l i=F r 

^if rs? =0, i = B r _i+2,..., J B r . (11) 

s=l 

respectively. As we demanded, the constraints in the Eq. 
© and Eq. (J7J are divided into independent constraints 
for each degenerate subspacc of p a . 

Also, from the construction of matrices X[ 1 we have 
K rs F r = 1/ Wm r . Using this relation along with Eq. (|lll) , 
we conclude that in Eq. © the term F r vanishes. Fi- 
nally, we can write 

2 2-. 

m n — 1 

ii P -n»ii 2 = ]r £4 

d 

+ J2{trC r C t r -trK r C r C t r K t r ), (12) 

r=l 

where C r is a (m 2 , — 1) x (?i 2 — 1) dimensional matrix 
with entries {c^ : i = B r —i + 2,...,_B r } and if r is a 
77i r x (m 2 — 1) dimensional matrix with entries {K rs i : 
s = l,...,m r ;i = B r -\ + 2, B r }. Moreover, notice 
that Eq. (fTU)) and Eq. (fTTJ) . in fact, are some constraints 
for each matrix K r independent of the others. So we 
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conclude that: 



max ||p- IP (p) || 



EE 4 

i=D j = l 



similar argument is true for an arbitrary m x n state p 
for which m r < 2. In the following, we illustrate this, in 
an example. 

Consider the 3 x n dimensional state p written in the 



Consequently, MIN can be written in the form 

d m 2 n 2 — 1 



i\r(p) = X>(p) + E E 



form of Eq. © as p = ^ ^ ® ^ + yj ± (x - §) (V3X 2 - 

+ ^max{trC r C t r -trK r C r C t r K t r ). X 3 +2X 5 + 2X?) <g> + <g> Vl Yi - =1 Ui X { <g> F x where 

g < x < and 2/1 and {t^i : £ = 1, ...,8} arc arbitrary 
real numbers chosen in a such a way that makes p to 
be a valid density operator. Also, in the computational 
bases, the operators {X; : i = 1, 8} are in the following 
form: X 1 = ^(|1)(1| - |2)(2|), X 2 = ^(|1)(1| + 12)(2| - 

2|3)(3|), X 3 = i(|l)(2| + |2)(1|), X 4 = 7 |(|1)(2| - 

|2)<1|), X 5 = -^(|1)( 3 | + |3)(1|), X, = 7 i(|l)(3|-|3)(l|), 

X 7 = ^(|2)(3| + |3)(2|) and X 8 = -§(|2)<3| - |3)(2|). 

Using Eq. we have N(p) < J2i=i^a- We now eval- 
uate our upper bound as follows. At first, we evaluate 
the reduced density matrix of p. The spectral decompo- 
sition of p a is p a = ei^2 s _ 1 \ki s )(ki s \ +e' 1 \ki)(ki\ where 
ex = 2x-l/3, e[ = -Ax + 5/3, \ku) = ^(|1) + 12) + 13)), 

|fci2> = ^(|1) - |2» and \h) = ^(|1) + |2) - 2|3)). We 
now construct the {X[ : i = 1, 9} which are used in the 
expansion of p in the form of Eq. (|5]). We construct {X^ : 
i = 1,...,4} as follows; X{ = ^(|fcn)(fcii| + |fc 12 )<fci 2 |), 

X 2 = ^(|fcn)(fci2| + |fci2>(fcn|), ^ = ^(|fcn>(fc 12 | - 
|*i 2 )(fcii|) and X' 4 = ^(|*ii)<*ii| - |fc 12 )(fc 12 |). In ad- 
dition we set X' 5 = \ki)(ki\. The operators {X- : i = 
6, ...,9} are constructed in the following form; X' 6 = 
^(|*ii>(*i| + |*i)(fc u |), X> 7 = ^(|fcii)(fci| - |fci)(fcu|), 



r=l i—D j=l 

N r (p) = trC r Cl - mm trK r C r C*K^.. (13) 

So, the evaluation of MIN for an m x n dimensional state 
p, is degraded to finding the MIN in the m r x n dimen- 
sional subspaces of state p. This method can reduce the 
calculations in the evaluation of MIN. 

As an important result, we obtain an upper bound for 
MIN from Eq. (TT31) . Following the same arguments in 
Ref. we can write a similar upper bound for each 
N r (p). So we have 



d rii r —m r 



m n —1 

N (p) ^ E E A - + E - 

i=D j=l 



(14) 



r— 1 i—1 



where {X r i : £ = 1, ...,m^ — 1} are the eigenvalues of the 
matrix C r C* listed in decreasing order. 

In the next section, we illustrate our previous results 
for 3 x n dimensional states p, as an example. 
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h(\ki2)(h\ + \ki){k 12 \), x> 9 = $(\k 13 )(k 



V2 



degenerate i. e. p a 
Thus, in this case, d 



3 x n DIMENSIONAL STATES 

Consider an arbitrary 3 x n dimensional p. As stated 
before, when p a is non-degenerate, the evaluation of N(p) 

is simple i. e. N(p) — Y1T=d Ej^i" 1 c ij- Now, we con- 
sider the case where one of eigenvalues of p a is doubly 

eiELil fc is)(M + ei|fci)(fci|. 
1, d' = 1, mi = 2 and D = 6. 

Using Eq. (USD we have 2V(p) = N^ + ^f^ 1 TLa c ij- 
Computing N\(p) is equivalent to computing MIN of a 
2 x n dimensional state p which its marginal p° is degener- 
ate. The analytical formula of MIN for an arbitrary 2 x n 
dimensional state p has been given in Eq. © of Ref Q . 
Using this equation we have Ni(p) = trC\C\ — \ m in 
where X m in is the smallest eigenvalue of matrix C\C\. In 
addition, in this case, the upper bound in Eq. (|14[) is 
equal to N(p). 

So, when one of eigenvalues of p a is doubly degenerate, 
our method leads to the exact value of the MIN. More- 
over, our upper bound in Eq. (|14[) coincide with the exact 
value of MIN which is, in general, better than the upper 
bound introduced in [3j, i. e. Eq. ([5])- Also note that a 



IfciXM)- 

Now, we write p in the form of Eq. © using {X[ : 
i = 1, 9}. We have p = (^X[ + ^X' b ) + 
Ei=i caX'i ® Yi where c n = ^t u + ^-t 31 - ±t 4 i + 
3^61 +^8i)i c 2 i = ^?(2i 2 i +tei —tgi), C31 = -^(—2^5i — 



tn + tgi), C41 



^5^31 + |*41 + 3(^6! + t 8 i), C51 



V2j 



v2 



[tei + t$i), cqi 



\/2 



V6 l 



V2 



(tn+tn), c 8 i = ^-(t 



i 6 i + tsi) and C91 = ^(*5i - *7l + *9i)- According to 
our previous descriptions (in the first paragraph of this 
section), we obtain Aq(p) = J2* =2 c fi- So AT(p) = — c\ x + 
y]i—2 c fi- As stated before, in this case our upper bound 
is equal to N(p). The difference between the upper bound 
introduced in Eq. ([5]) and our upper bound is equal to 
g(^*3i + V2t 61 + V2t 81 ~ ^4i) 2 - So our upper bound 
is better than the upper bound introduced in Eq. (|5|). 

Now consider the case that p a is fully degenerate i. e. 
p a have three identical eigenvalues. In this case, we have 
d = 1, d' = 0, mi = m = 3 and so Eqs. flO]). fTTj). 
(TT3]) are the same as Eqs. ©, ([7]), (O, respectively. We 
use the notation of Eqs. ([5]), (|7I3|) for simplicity. When 
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p a is fully degenerate, then p a = — . So the constraint 
J2k^kP a ^t = /° a sct n0 conditions on {II^}. Therefore, 
in the optimization problem of Eq. only the con- 

straints in the Eqs. (|6|), (O are present. Using the {X{\ 
represented in the previous paragraph, we can obtain the 
constraints on {aki} which can be derived from Eqs. (JSJ), 
([?])• These constraints are represented in the appendix. 

In order to evaluate MIN, we must find the minimum 
of tr ATT 1 A 1 with the constraints obtained from the Eq. 
© and Eq. ((7|) in the previous paragraphs. But, if we 
put additional constraints on the {aki} simplifying the 
previous represented constraints, then we can find some 
analytical lower bounds for MIN. For example if we put 
in = 0,21 = 0, a%2 = a.21 = and a 2 6 = o 28 = then we 
obtain a 3 i = a 32 = 0, a 25 = a 27 = a 13 = a 23 = a 33 = 

~ 2 3^ i a 15 = a 17 = 0,35 = a 37 = &16 = a 18 = T^g, 

a 36 = a 38 = an d ai4 = a 2 4 = a 3 4 = and thereby 
we obtain a lower bound for MIN. For every state p with 
{tij = t2j = tij = : j = l,...,n — 1}, this lower 
bound is equal to the exact value of MIN. So depending 
on the coefficients tij in the expansion of the the state 
p, this lower bound can be a good lower bound for MIN. 
According to this method, we can find good lower bounds 
for MIN of the other classes of states. 



CONCLUSIONS 

We have studied the evaluation of MIN for an arbi- 
trary m x n dimensional bipartite density matrix p. We 
have shown that according to the degeneracy of p a , ex- 
panding p in the suitable basis, the evaluation of MIN 
for an m x n dimensional state p, is degraded to finding 
the MIN in the rrii x n (rrii < m) dimensional subspaces 
of state p. This method can reduce the calculations in 
the evaluation of MIN. Moreover, for an arbitrary m x n 
state p for which rrii < 2, our method leads to the exact 
value of the MIN. Also, we obtain an upper bound for 
MIN which can improve the ones introduced in Ref. 
In addition, we have explained the evaluation of MIN 
for 3 x n dimensional states and introduced some lower 
bounds for MIN when p a is fully degenerate. 
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APPENDIX 

In this appendix, we represent the constraints on {aki} 
which can be derived from Eqs. ©, ([7]). From Eq. 
for k' = k = 1, 2, we obtain (we recall that aio = a 2 o = 

75) 

T 8 a 2 - 2 

Z^i=l u ki ~ 3 

\ a ki = ^Ok\Ok2 + ^75 {a 2 k5 + a| 6 - a| 7 - a| 8 ), 

\a k 2 = 7g(afel- a L+ a fc3 + a L-5( a fe5+ a fe6 + a fe7 + a fc8))' 

\okz = "-^afc 3 afe 2 + --^(a^a^ + a^aks), 

\om = -^OkAOk2 + -^{ — Ok50-m + Ok60k7), 

|afc5 = -J=0k50k2 + -^{oklOkh + OkzOkJ — a^Ofcs), 

\ake = ^Ok&Ok2 + -^{akiOke + afc30fc8 + a^a^), 
\ok7 = 7^=a/c7afc 2 + t^(— Okiaki + afe 3 afc5 + a^a*^), 

|afe8 = -^Ok80k2 + T v 7 2 = ( — a fcl fl fc8 + afe30fc6 — OkA,Okb)- 

In addition, if we put k = 1 and k' = 2 in the Eq. (|6|) , 
then the following constraints on {aki} can be obtained. 

V^8 _ -1 

2_,i=l a li a 2i — 

i(an + a 2i ) + -^(ana 22 + a i2 a 2 i) + 5-75(015025 + 
awa 2 6 - ai 7 a 27 - ai 8 a 28 ) + 5775(20140^ - 2ai 3 a 24 + 
016025 - Oi 5 a 26 + oi 7 a 28 - oi 8 a 27 ) = 0, 

|(ai 2 + a 22 ) + 7^(ana 2 i - ai 2 a 22 + a i3 a 23 + ai 4 a 24 - 

50l5 a 25 — ^016026 — 5^17027 — \a\%02%) + 7^7g(ai6 a 25 — 
ai5026 + &18027 - Oi 7 a 28 ) = 0, 

i(a i3 + a 23 ) + --)g(ai 2 a 23 + a i3 a 22 ) + 5775(0150^7 + 

ai7»25 + ai6 a 28 + a 18 a 26) + 775 (on0 2 4 — 014^21 — 5ai5 a 28 + 

^ai8a 2 5 + ^ai6a 2 7 — ^017026) = 0, 

\{a u + a 24 ) + 7^(ai 2 a 24 + ai 4 a 22 ) + ^(-ais^s - 
ai8d25 +ai6«27 + ai 7 a 26 ) + (ai 3 a 2i - ana 23 - iai 5 a 27 + 
^ai7 a 25 + 2 1 ai8 a 26 — ^ai6a28) = 0, 

i(ai 5 + 025) - 5^(012025 + 015022) + 7^75(011025 + 
ai5O 2 i+ai 3 a 2 7 + ai7a 23 -ai40 2 8-ai 8 a 2 4) + 5-j 7 5(aiia 26 - 

ai6fl2i + 013028 - ai8<223 + ai4a 27 - ai 7 a 24 + V3ai 2 a 26 - 
V / 3ai6a 22 ) = 0, 

i(ai 6 + a 26 ) - 5-j7g(ai2a26 + 016022) + 5775(011026 + 
ai60 2 i+ai 3 a 28 +ai 8 a 23 +ai4a 2 7+ai 7 a 24 ) + 77^(-aiia 2 5 + 

ai5a 2 i - ai 3 a 27 + ai 7 a 23 + ai 4 a 28 - ai 8 a 24 + \/3ai 5 a 2 2 - 
v / 3ai 2 a 25 ) = 0, 

|(ai7 + a 27 ) - 57yg(ai 2 a 2 7 + ai 7 a 22 ) + 5-y5(-aiia 27 - 
017021 +ai3025+oi5a 23 +ai4a 2 6+ai6a 2 4)+5-y5(-oiia 28 -|- 
oi8a 2 i + ai 3 a 26 — ai 6 a 23 — ai4a 25 + ai 5 a 2 4 + -y/Sa^a^ — 
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\/3ai8 a 22) = and 

i(ai 8 + a 2 s) - j^(ai2a 2 s + a w a 2 2) + ^(— 011028 ~ 
018^21 +ai 3 a 26 + ai 6 a 2 3 - a u a 2b - ai 5 a 2i ) 4-^(a n a 2 7 - 

ai7«2i - ai 3 a 2 5 + ai 5 a 2 3 - ai 4 a 26 + a w a 2 A + V3ai 7 a 22 - 
v / 3ai2a27) = 0. 

From Eq. ([7]), for i = 1, m 2 , wc obtain 0,3^ = — (ai»+ 
a 2 i). Using this relation, it is concluded that by choosing 
k, k' = 3 in the Eq. ©, no additional constraint can be 
obtained. 
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